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One-loop mass shifts to the classical masses of stable kinks arising in a massive non-linear S 2 - 
sigma model are computed. Ultraviolet divergences are controlled using the heat kernel/zeta function 
regularization method. A comparison between the results achieved from exact and high-temperature 
asymptotic heat traces is analyzed in depth. 
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O ■ I- INTRODUCTION 

bJQ 1 In a seminal paper, Olive and Witten [l[ linked extended supersymmetric theories to BPS solitons by showing that 
the classical mass of these stable lumps agreed exactly with the central charge of the extended SUSY algebra. The 
subsequent issue concerning BPS saturation at one-loop (rather than tree) level has proved to be extremely subtle, 
prompting a remarkable amount of work over the last twelve years. See, e.g., and References quoted therein to 
find an in-depth report on these developments. 

A new actor entered the stage when in Q a Stony Brook/Wien group computed the one- loop mass shift of the 
supersymmetric CP 1 -kink in a N = (2, 2) non-linear sigma model with twisted mass. Kinks of several types in 
massive non-linear sigma models were, however, discovered earlier, see [H, Q. In Reference H, three of us 

found several families of non-topological kinks in another non-linear sigma model: we chose § 2 as the target space and 
considered the case when the masses of the pseudo-Nambu-Goldstone particles were different. The 0>(2)-symmctry of 
(— i ' the equal-mass case is explicitly broken to Z2 x Z2 and the SO(2)-families of topological kinks of the former system are 
l — 1 , deformed to the four families of non-topological kinks arising in the second system. The boundary of the moduli space 
of non-topological kinks in the last model is formed by a pair of topological kinks of different energy. The analysis 
of kink stability in the massive non-linear § 2 -sigma model performed in [9( allowed us to calculate the one-loop mass 
shifts for the topological kinks by using the Cahill-Comtet-Glaubcr formula These authors showed that the 

ly-j I one-loop mass shift for static solitons can be read from the eigenvalues of the bound states of the kink second-order 
£NJ ■ fluctuation operator and the threshold to the continuous spectrum when this operator is a transparent Schrodinger 
operator of the Posch- Teller type. This is the case of the topological kinks of the massive non-linear § 2 -sigma model 
\q when a parallel frame to the kink orbits is chosen to refer to the fluctuations. 

The aim of this paper is to offer another route for computing the one-loop kink mass shift in order to unveil some of 
0^ . the intricacies hidden in this subtle problem. We shall follow the method developed in References [l2j and [l3| based 
■ on heat kernel/zeta function regularization of ultraviolet divergences. See also the lectures [l4|, where full details 
can be found. Because the spectrum of small kink fluctuations in our system can be identified analytically, we are 
able to give the exact answer for the mass shifts. We shall also show, however, how to reach approximately the same 
result using the coefficients of the heat kernel asy mptotic expansion. The interest of this calculation is that a formula 
belonging to the class of formulas shown in [17[ will be derived. The importance of this type of formula lies in the 
fact that it can be applied to obtain the one-loop mass shifts of topological defects even when the spectrum of the 
second-order fluctuation operator is not known; for instance, in the case of two-component topological kinks: see [12], 
[l3l ]. Similar formulas work even for Abelian gauge theories in (2+l)-dimensions and thus the mass shifts of self-dual 
Nielsen-Olesen vortices and semi-local strings can be calculated approximately, see [l8[, [Tj|, and poj . 

To end this brief Introduction we simply mention that interesting calculations have recently appeared addressing 
one-loop kink mass corrections and kink melting at finite temperatures in the sine-Gordon, CP 1 , and A^ 4 models in 
a purely bosonic setting, see [l5j . 

The organization of the paper is as follows: In Section § .II, we introduce the model and explain our conventions. 
In Section § .III, the perturbative sector as well as the mass renormalization procedure are discussed. Section §.IV is 
devoted to the analysis of the stable topological kinks in this system. The second-order kink fluctuation operator is 
obtained, placing special emphasis on its geometric properties. In Section §.V, the one-loop mass shift is computed 
using the heat kernel/zeta function regularization method. Section §.VI offers a comparison of the exact result obtained 
in §.V with the approximation reached from the high-temperature asymptotic expansion. Finally, a summary and 
outlook arc offered whereas two Appendices containing some technical material are included. 
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II. THE (l+l)-DIMENSIONAL MASSIVE NON-LINEAR § 2 -SIGMA MODEL 

The action governing the dynamics of the non-linear § 2 -sigma model and the constraint on the scalar fields are: 

S[^,M = fdtdxUg^j^^} > ■ (1) 

The scalar fields are thus maps, <j> a (t, x) G Maps(R 1-1 , § 2 ), a = 1, 2, 3, from the (1 + l)-dimensional Minkowski space- 
time to a § 2 -sphere of radius R, which is the target manifold of this non-linear sigma model. Our conventions for R 1,1 
are as follows: x^ G R 1 ' 1 , fi = 0, 1 with x° = t,x x = x and g 1 *" = diag(l, —1). Then x^ ■ x^ = g^ x^Xy = t 2 — x 2 and 

d f 8 \_ „„ d 2 d^ _j¥_ 

dxn \dx^ J ^ dx^dx u dt 2 dx 2 

The infrared asymptotics forbids massless particles in (1 + l)-dimensional scalar field theories, s ee fig . We therefore 
include the simplest potential energy density that would be generated by quantum fluctuations [2l|: 

F(01, 2 , 3 ) = - (a\ 2 + a\ 2 + a\ $) 



1. Solving 3 in favor of 0i and 02, sg(03)03 = \JR 2 — <f>( — 0|, we find 

G 1 \ d <f>l d 4>l , d(f> 2 902 , (01^01+02^02) (01^01+02^02) t * r 



where 



F S2 (0i, 02 ) = i ((a 2 - a 2 ) 2 + (a 2 - a 2 ) 2 + const.) ~ y 2 (t, s) + y 2 ,(i, s) 



with A 2 = (a 2 — a 2 ), -f 2 = (a 2 — a 2 ), X 2 > -f 2 . The masses of the pseudo-Nambu-Goldstonc bosons are 
respectively A and 7. 

2. Interactions, however, come from the geometry: 

(01<9 M 01 + 02l9 M 2 ) (01<9 M 01 + 02<9 M 02) 



v/i? 2 - 2 - 2 v/i? 2 - 2 - 2 

and -j^z is a non-dimensional coupling constant. 

In the unit natural system, h = c = 1, the dimensions of fields, masses and coupling constants are respectively: 
[0 a ] = 1 = [/?], [7] = M = [A]. We define non dimensional space-time coordinates and masses: 



rM ^ 2 _ gj " gj _ 7 2 

A ' a 2 -a 2 A 2 

to write the action and the energy in terms of them: 



< a 2 < 1 



1 f MA J ^ 9 fo J. ^ d ^ _1_ Wj+W) (01^01+ 02^02) , 2 .. , 2,2,, x I , s 

dtda; { — ■ — + — ■ — + — --^fosj-g ^(t.arU (2) 



^-01 



E—-fd ' ^ < f d< ^ 2 \ 2 i (01^01 + 02l9 t 2 ) 2 , fd(j)i\ 2 i fd(j} 2 \ 2 , (01 <9 X 01 + 02 ^02) 2 



2 7 ""~ Ha/ ' \ dt J R 2 -cj) 2 -<p 2 \ dx J \ dx J R 2 - <f>l - 0| 

+ 2 (t,x)+a 2 -0 2 (t,a;)} 

There are two homogeneous minima of the action or vacua of our model: 0f = 2 / =0,0^ = ±i?, the North 
and South Poles. Choice of one of the poles to quantize the system spontaneously breaks the Z2 x Z2 x Z2 symmetry 
of the action ©, a -> (-l) <5 » b o , a, 6 = 1,2,3, to: Z 2 x Z 2 , Q -» (-1) 5 °"0^, a, = 1,2. Therefore, the 
configuration space C = {Maps(IR, S 2 )/^ < +00} is the union of four disconnected sectors C = Cnn U ^ss U ^-ns U ^-sn 
labeled by the vacua reached by each configuration at the two disconnected components of the boundary of the real 
line: x = ±00. 
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III. MASS RENORMALIZATION 



The field equations 



□ c 



R 2 - ELl §a4>a 

v — > 2 
>2 La=l <? 



0=1 



R 2 - Ea=l 

become linear for small fluctuations, G a (x fl ) = 
n5Gi(t,x) + SGi.{t,x) = 0(5G a 6Gp) 



ELi ^ad^a J2l=i ^pd^cjip 







(JGa^a;^), around the vacuum: 

aSG 2 {t,x)+a 2 6G 2 {t 1 x) = 0(SG a SG p ) 



(3) 



We shall need the Feynman rules only for the four-valent vertices. Besides the two propagators for the (pseudo) 
Nambu-Goldstonc bosons there are three vertices with four external legs. The derivatives appearing in the interactions 
induce dependence on the momenta in the weights. This also affects the sign and the combinatorial factors. Naturally, 
there are many more vertices in this model, but we list only the vertices that contribute to the self-energy of the 
Nambu-Goldstone bosons up to one-loop order. 

Table I: Propagators 



Particle 



Field Propagator 



Diagram 



Nambu-Goldstone Gi(a;' i ) 



k 2 - k 2 ~ 1 + ie 



Nambu-Goldstone Giix 11 



fcp — k 2 — a 2 + ie •-- 



Table II: Fourth-order vertices 



Vertex Weight Vertex Weight Vertex Weight 




A. Plane waves and vacuum energy 

The general solution of the linearized field equations ([3]) governing the small fluctuations of the Nambu-Goldstone 
fields is: 



SG 1 (x ,x) = \ ■ ^ 2 l i(k) {ai{k)e- ik ° x ° +ikx + a{(k)e ik °*°- ikx } 

5G 2 (x ,x) = \ ■ J)Y, J n {^(l)e i90X °- iqx +^(g)e-^°+^} 
* V I y/ 2aj 2 {q) 
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where ko = u>i(k) — \Jk 2 + 1, qo = 012(0) = V / q 2 + & 2 , and the dispersion relations k 2 — k 2 — 1 = 0, q 2 — q 2 
hold: 




f #011 





( 


V 


^022 y 





-2- 4- 1 





We have chosen a normalization interval of non-dimensional "length" I = XL, I = [—4, 5], and we impose PBC on 
the plane waves so that: k ■ I = 2xni , q ■ I = 2im 2 with nx,U2 € Z. Thus. A'o acts on L 2 = © a=1 L 2 (S 1 ), and its 

1 o\ 



spectral density at the I — > 00 limit is: pk (fc) = ^ dna 

' " "ST 



2tt 



1 



From the classical free (quadratic) Hamiltonian 



dx 



( dSGi d6Gi d8(h dSGA ( d8G 2 dSG 2 d8G 2 dSG 2 



\ dxo dxo 



dx 



dx 
A 



\ dxo dxo 



dx 



dx 



+ Sd ■ 5G 1 + a 2 5G 2 ■ 5G 2 } = £ £ - [cj a (k)(a* a (k)a a (k) + a a (k)a* a (k))} 

k a=l 

one passes via canonical quantization to the quantum free Hamiltonian: 

[a a (k),al(q)} = S a0 5 kq , = £ A wi(fe) ^(fc)ai(fc) + + w 2 (fc) ^(fc)a 2 (fc) + ij 

The vacuum energy is: 

a a (fc)|0 ; y} = 0, Vfc,V« , AS = (0;^|^ 2) |0;y) = ^ Wl (fc) + ^c l ; 2 (A : ) = ^Tr L2 A| 



B. One-loop mass renormalization counter-terms 

There are four ultraviolet divergent graphs in one-loop order of the ^.-expansion contributing to the Gi(x^) and 
G 2 (x fl ) Nambu-Goldstone bosons self-energies: 



• Self-energy of G 2 



p p p p 

12 1 



2i f d 2 k iiPn + k^p* 1 2i f d 2 k i( Pfl + fc M )p" 



R 2 J (2tt) 2 (p M + kjk" - 1 + is R 2 J (2tt) 2 (p^ + k^ - a 2 + ie 
2i f d 2 k i 2i f d 2 k i 



R 2 J (2ir) 2 k^-l + ie R 2 J (2tt) 2 fc M fc^ - a 2 + is 

• Self-energy of G\ 

P+k/^~~\ p+k, 
2ia 2 _ 2io 2 1^1 



1(1) + _.J(0= v 



i? 2 w i? 2 v ; | — 

P P 
2icr 2 r dk 1 2ia 2 f dk 1 




i? 2 y 4^ VFTT ^ 2 y 4^ Vfc 2 + 
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where we have computed the kg integrations by using the residue theorem. We only show this step explicitly 
in the computation of the self-energy of SGi because it suffices to point out how to regularize these divergent 
integrals by means of spectral zeta functions. The rcgularization just mentioned will be performed later in 
Section §.V D. 

The p^p^ factor becomes constant when the momentum is put "on shell" in the external legs, p^p^ = 1, 
p^p^ = u . This process gives us the mass renormalization counter-terms. The Lagrangian density of counter- 
terms shown in Table III must be added to cancel the above divergences exactly. We also show the vertices 
generated at one-loop level. 



Table III: One-loop counter-terms 



Diagram Weight 

y -f /(- 2 )) 11 



IV. ISOTHERMAL COORDINATES AND TOPOLOGICAL KINKS 

In this Section we shall use the isothermal coordinates in the chart S 2 — {(0,0, —R)} obtained via stereographic 
projection from the South Pole: 

1 = 01 = gfr 2 = 02 = R4>2 (4) 

1 + f i? + sg(0 3 )^ 2 -0i-02 ' 1 + f R + sgfoW 

The geometric data of the sphere in this coordinate system are: 

j 2 4i? 4 f 1 2\ I 1 2\ 4 i? 4 

dg = (jR 2 +x l x l +x 2 x 2)2 . 9ll(X ,X )= 922(X ,X ) = 

r\iix\x 2 ) = -r^ 2 (x 1 ,x 2 ) = rf 2 (x 1 ,x 2 ) = r| 1 (x 1 ,x 2 ) 



(R 2 + x'x 1 + x 2 x 2 ) 2 y VA ^ ' > (R? + x V + x 2 x 2 ) 2 

2X 1 



R 2 + x 1 x 1 +x 2 x 

rUx\x 2 ) = -r 2 11 (x\x 2 ) = rl 2 ( x \ x 2 ) = rl 1 ( x \ x 2 ) = 



2 



R\22(x\x 2 ) = -Rli2(x 1 ,x 2 ) = -Rl2i(x 1 ,x 2 )=R%u(x\x 2 ) 



R 2 + x^ x +x 2 x 2 

-AR 4 



(R 2 + x 1 x 1 +x 2 x 2 ) 2 

The kinetic and potential energy densities read: 

nx\x 2 ) = ?f M 2 2 , 2 • {dtx'dtx 1 + d tX 2 d tX 2 ) 

(R z + X X + X X J 

V{ X \ X 2 ) = a 2 , 2 • (d.x^x 1 + d xX 2 xX 2 + X 1 X 1 +a 2 X 2 X 2 ) 

(R z + x X + X X J 

From the action S = J d 2 x [T — V] one derives the field equations: 

□X' + rj^xW + ^X 1 + o 2 5l X 2 - 2{5\ X 1 + S 2X 2 ) fff^ 2 = . ( 5 ) 

^ i X X i X X 

which for static configurations reduce to: 

d 2 X l -™ dx> X k ! 2 2 i 1 , « 2\ xV + cr 2 x 2 X 2 n ,^ 

- -r-2- -T Jk — — +S\x +cr S 2X -2(6\x + S 2X ) = . (6) 

ax- 2 J ax ax R z + XX + X X 
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A. Topological K kinks 

We try the x 1 = orbit in ([6]) and reduce this ODE system to the single ODE: 

d 2 X 2 2X 2 d X 2 d X 2 2 2 f X 2 X 2 \ 

~ a X 1 - 2 n9 ■ „ , . (7) 



dx 2 R 2 + x 2 X 2 dx dx \ R 2 + x 2 X 2 

Xk(x) = , X 2 K{x)=±Re ±a{x - Xo) , (8) 

are solutions of ([7]) of finite energy: 

/OO /? 2 fT 2 
d.T = 2XR 2 a . (9) 
- oo cosh (a(x - x j) 

In ([8]), is an integration constant that sets the kink center. The kink field components in the original coordinates 

0f(x)=O , 0f (x) = f , <f%(x)=±R\wib{o(x-x Q )] 

cosh[cr(a; — xoj\ 

are either kink-shaped, or bell-shaped, It is clear that the four solutions (JSJ) belong to the topological sectors 
Cns or Csn of the configuration space. Lorentz invariance tells us that 

Xk(*)=0 , XkW = ±iZ exp[±a( *Z^L - x )] (10) 

V 1 — u 

are solitary wave solutions of the full field equations ((5]). 

B. Second-order fluctuation operator 

Let us consider small kink fluctuations: 

X(x) = Xk(x) +i](x) , r)(x) = (i] 1 (x),r] 2 (x)) 

Here, xk{x) = (xk i x )i Xk( x )) is the kink solution and r](x) = ^(x)-^- + r/ 2 (x)-^ e T(T§ 2 ) are vector fields along 
the kink orbit - expressed in the orthonormal basis A-) = <5 U of T§ 2 - giving the small fluctuations on the kink. 



From the tangent vector field to the orbit x'k( x ) = ~SFd§T + "S^af 7 ' ^ ne covariant derivative, and the curvature 



v d x % ' 9 X 3 

he tangent vector field to the orbit Xk( x ) ~ 

tensor 

V x ' K v(x) = (r}' i {x)+T i jk { X K)rf{x)x'k{ x )) ^ R(x'k,v)x'k = X%( X W ' { x )x' k { x )R\, k {x K )^ 

we obtain the geodesic deviation operator: 

5W = ^^*) > ^(x)+R{Xk,v)x'k(x) ■ 
We also need the Hessian of the "mechanical" potential U(x 1 , X 2 ) = — ^(x 1 , X 2 ) 

( d 2 U dU \ 3 

V,gradL/(x) = [^-{ XK ) - T% (xk)^ ( X k) J <7^ 

The second-order fluctuation operator around the kink xk is: 

-£ 2 77 



A(iT)r/(x) = - 



dx 2 



(x) + R(x' K , v)x'k + V„grad[/ (x) 



(11) 
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C. Small fluctuations on K kinks 

Application to the K kink Xk{x) = (Xk( x ) = ®>Xk( x ) = Re~ ax ) gives: 
A(K)t] = 



1 2<t(1 - tanh<7x))-j- - (l - 2<r 2 + 2(T 2 tanhcra:) ?/ ) — 
1 2ct(1 - tanhax))— — + a 2 (1 - 2tanh<7x) rf ' 



c?x 2 " dx y ' 1 J dx 1 

(12) 



dx 2 dx J d\ 2 

The second-order fluctuation operator in the orthonormal frame is a second-order differential operator that has first- 
order derivatives both in the direction of the kink orbit, -^s, and the orthogonal direction to the orbit f^r- 

Alternatively, we can use a parallel frame, fi(x) = fj}(x)-^r + n 2 (x)-^, along the K kink orbit: 



% +r i jk (xK)x'^ k 



&r+<j(l- tanh)^ 1 (x) = => [i 1 (x) = l + e~ 2ax 
^ +<t(1 -tanh)^ 2 (a;) = => [i 2 (x) = 1 + er 2ax 



In this parallel frame the vectors of the basis ^ (x)-^ point in the same directions as but their moduli vary along 
the kink orbit: 



Writing the fluctuations in this frame, rj(x) = (x) /i 1 (x) g^- + ^ 2 (x)/i 2 (x) we find 



kit* us( 2<j2 s+\\ 9 i, J d2 ¥ ,i 



2d 2 x A 9 



cosh 2 crx / dx 2 



)r £o ■ (13) 



In the parallel frame the second-order fluctuation operator is a transparent (reflection coefficient equal to zero) Posch- 
Tcllcr Schrodingcr operator both in the parallel and orthogonal directions to the kink orbit. 

This analysis is deceptively simple: acting respectively on r) l (x) = (l + e~ 2<TX )t; 1 (x) and r] 2 (x) = (l + e~ 2(TX )t; 2 (x) the 
terms with first-order derivatives in (|12|) disappear and (1 + e~ 2ax ) factors out, leaving very well known Schrodingcr 
operators acting respectively on ^(x) and £, 2 (x). The key point is that the differential operators in (fT2|) and (fl3|) 
share the eigenvalues although their eigenfunctions differ by the /i l (x) factors. The spectral functions associated are 
thus identical and it seems wise to use the best known form. What we have shown here is the geometrical meaning 
of the n % {x) factors: they provide a parallel frame along the kink orbit. 

D. The spectrum of small kink fluctuations 

Changing from vector to matrix notation, 

d ( l\ _d_ 

dx 1 ^ \ J ' d X 2 



we now use the differential operators of formula (|13|) to write the linearized field equations satisfied by the small kink 
fluctuations in the parallel frame: 

X 1 (t,x)= X 1 K (x)+fi 1 {x)SK 1 (t,x) , X 2 (t,x)=xl(x)+LJ 2 {x)8K 2 {t,x) 
^ + M=fl , ^+K 22 SK 2 = . 

Therefore, the eigenfunctions of the differential operator 

d 2 , I 2a 2 n \ 





° 


-( 


V o 


k 22 j 





K =\ n ^ ) = ( &X o COSh - _d 2 + (J 2_ 2 a 2 (14) 

dx 2 cosh^ ax / 



provide the general solution of the linearized equations via the separation ansatz: 8Ki(t, x) = .9i(i)C 1 ( a; ); 8K 2 {t, x) 
cj2(t)t; 2 (x). The eigenvalues and eigenfunctions of K are shown in the following Table: 
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Eigenvalues 


Eigenfunctions 


Eigenvalues 


Eigenfunctions 


£ i-*z = 1 - a 
e\{q) = a 2 q 2 + 1 


V J V / 

( fq{x) = e^™ (t&nhax - ik) \ 

V J 


e = 

ei(Jfe) = a 2 (k 2 + l) 


( o W \ 

( ) 

V /fcW = e lfc,ra (tanher:c - ik) J 



The spectrum of K22 contains a bound state of zero eigenvalue -the translational mode- and a branch of the 
continuous spectrum, with the threshold at e|(0) = a 2 . Speci^n also is formed by a bound state of positive eigenvalue 
and a branch of the continuous spectrum starting at £i(0) = 1. Periodic boundary conditions in the [— |, |] interval 
require: 

aq ■ I + S\ (q) = 27m 1 , ok ■ I + 82 (k) = 27m 2 , 711,712 6 Z , 
such that the phase shifts and the induced spectral densities are: 

Si(q) — 2arctan- = S(q) , ^(fc) = 2arctan— = S(k) 
q k 



PK ^ q) = h{ al+ i^ [q) ) ' PK ^ {k) = h{ al+ it {k) ) ■ (15) 

In sum, K also acts in the Hilbert space L 2 = @ 2 a=1 L^S 1 ), and its spectral density in the limit of very large radius 
of the circle is: 

-«-(T;H(-s»)(ij) ■ 



V. ONE-LOOP SHIFT TO THE CLASSICAL K KINK MASSES IN THE MASSIVE NON-LINEAR 

S 2 -SIGMA MODEL 



A. Zero-point kink energy 



The general solution of the linearized field equations governing the small kink fluctuations is: 



6Ki(x 0l x) 



A x . 



O i\/l — a' 2 xo 



— i^/l — a' 2 xo 



fl-a 2 (x) 



SK 2 (x 0l x) 



1 1 

2 V2Vl - o 2 

\ ■ J]J2 {A 1 {k)e- i ^ k ^fl{x) + Al(k)e^W^f^(x)} 

\ ■ y=^y {A2(k)e-^ k ^f 2 {x ) + A* 2 {k)e^ k ^ft 2 {x)) 



Note that the zero mode is not included because only contribute to quantum corrections at two-loop order. In the 
orthogonal complement to the kernel of K22 in (Ba=i -£'a(^ 1 )i the eigenfunctions of the K operator satisfying PBC 
form a complete orthonormal system. Therefore, the classical free Hamiltonian 



H< 2 > = E 



A ^ f dSKg dSK a 
V dxo dxo 



SK a ■ K aa ■ SK Q 



Vl-* 2 {Al^A^ + A 1 _ a ,A* 1 _, 2 ) {e a (k)(A* a (k)A a (k) + A a (k)A* a (k))} 



k a=l 
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can be written in terms of the normal modes of the system in the quadratic approximation. From this expression, one 
passes via canonical quantization, [A a (fc), j4t(g)] = 5 a f}8k q , [A 1 _ a .2, J Aj_ o . 2 ] = 1, to the quantum free Hamiltonian: 

= e+ \Vi-<j 2 (iUUi-^ + 1) + \J2 E Ui(k)A a (k) + 1) 

\ ' k a=l L V / 

The kink ground state is a coherent state annihilated by all the destruction operators: 

A a (k)\0;K) =A 1 _ (7 2\0;K) =0 ,Vfc, Vq , 4> a (t, x)\0; K) = <fi (x)\0; K) 
The kink ground state energy is: 



A + A£ = (0; A|iJ (2) |0; A) = 2XR 2 a + ^1 - a 2 + £a ^ = 2AA 2 cr + ^Tr^ A3 



(16) 



k a=l 



B. Zeta function regularization and Casimir kink energy 

I i 

Both Tt^Kq and Tr^K a are ultraviolet divergent quantities: one sums over an infinite number of eigenvalues, 
and a regularization/renormalization procedure must be implemented to make sense of these formal expressions. We 
renormalize the zero-point kink energy by subtracting from it the vacuum energy to define the kink Casimir energy: 



AE C = AE - AE Q = - 
2 



The subtraction of these two divergent quantities is regularized by using the associated generalized zeta functions, 
i.e., we temporarily assign to AE C the finite value: 

AEC{S) = 2 (aJ) S [TrL2A ^ S " TrL2Ar] = 2 (0 [CA ' (S) ~ C * (S)] 
at a regular point of both (k(s) and (k ( s )- Here, 

Ck(s)= A_S ' Cx (s)= E A o S ; A e SpccA , A G SpccA , seC 

Specif SpccATo 

are the spectral zeta functions of A and Kq, which are meromorphic functions of the complex variable s. An auxiliary 
parameter fi with dimensions of inverse length is used to keep the physical dimension right and we shall go to the 
physical limit AE C = lim 8 _,._i AE c (s) at the end of the process. 



C. Partition and generalized zeta functions 



Because analytical information about the spectrum of A is only available at the the limit of large I (bound state 
energies, phase shifts and spectral densities) it is better to consider first the partition or heat functions: 



Tr 



L 2 e 



-0K O 



2^ 



I 



(e- 



(3e 



Note that here we have replaced fc and q defined in Section §3.1 by ak and erg for a better comparison between the 
spectra of A and A. The PBC spectral density of A is thus obtained by replacing A by 7. The A-heat function is 
also expressed in terms of integrals over the continuous spectrum at the I = 00 limit, rather than infinite sums. The 
integrals, however, must be weighted with the PBC spectral densities: 
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The error and complementary error functions of j3 = jjAf , a fictitious inverse temperature, arise and the asterisk 
means that we have not included the zero mode because zero modes do not enter the one-loop formula (|16p . 
The generalized zeta functions are Mellin transforms of the heat functions: 
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We indeed find meromorphic functions of s with poles and residues determined from the poles and residues of Euler 
r(s) and Gauss hypergeometric 2 Fi [a, b, c, z] functions [22| . 

In the APPENDIX I we show that the kink Casimir energy in the physical limit s = — \ is the divergent quantity: 



. „n Xa 

AE C = 

2tt 



16 



lim 1 + 2 In & + in - 2 + a ff [1, 0, ^ 



1 



(17) 



where 2-pf ' 1 ' [g, 0; f, ~ ^ 2 ] is the derivative of the Gauss hypergeometric function with respect to the second 
argument. 



D. Zeta function regularization of the self-energy graphs and kink mass renormalization 



It remains to take the effect of mass renormalization into account. The contribution to the kink energy due to the 
mass renormalization counter-terms is: 



AE MR = _ x 



dxCc.T.^f 1 ,^ 1 ) = A^[/(l) +/(^ 2 )] / dx^(x)^(x) = 2Xa[I{\) + I(a 2 )] 



In the normalization interval of length I the integrals become infinite sums 
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that can be regularized by using zeta functions: 

m = -\ km^r^^^c^Gs + i) 

\xL s -»-i X z T(s) 



CK022 (S + 1) 



such that [22 
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1 
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In the APPENDIX I it is proved that the physical limit s = -\ is also a pole of AE MR (s): 
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Xa 
2^ 
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(18) 



The divergent terms in AE C (fl7|) and AE MR ([22]) . as well as the ^-dependent terms, cancel each other exactly and 
the one-loop K kink mass shift is: 
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(19) 



In formula ()19() we have also written the result found in our derivation a la Cahill-Comtct-Glauber of the quantum 
correction, see [|[. The heat kernel/zeta function result is ——/(c) whereas the CCH formula leads to — —g(a), 
where 



/(*) 



1 + 7T 2F i 



(0,1,0,0) 



bo,-,- 



5(c) 



Vl-c 2 



Despite appearances, /(c) and g(c) are identical functions of c 6 [0,1], as the Mathematica plots in the Figure 1 
show. This is remarkable: there is no mention about the analytic identity between the functions /(c) and <?(c) in the 
ample Literature on special functions. Nevertheless, they trace identical curves as functions of c. 
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Figure 1: Graphics of f(u) (left), g(o~) (center), and f(a) and g(a) plotted together (right). In the Figure, a is labeled as s in 
the abscissa axis. 



VI. HIGH-TEMPERATURE ASYMPTOTIC EXPANSION 



The exact heat or partition function can be written in the form: 
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tr 
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Tr 
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where "tr" means trace in the matrix sense. There is an alternative way of computing this quantity by means of a 
high-temperature asymptotic expansion. Although we have the exact formula in our system, we shall also perform 
the approximate calculation, which is the only one possible in other systems in order to gain control of this second 
approach in this favorable case. 

In the APPENDIX II it is shown how the coefficients of the power expansion of the K-he&t trace 



Tre 



-0K _ 



1 
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(20) 



the Seeley coefficients c n (K), are obtained through integration of the Seeley densities over the whole line. The densities 
satisfy recurrence relations tantamount to the heat kernel equation starting from a general potential U (x) . In our 
problem we must solve the recurrence relations between these densities for the potential U(x) — — — ^ ax , essentially 
the same potential as for the sine-Gordon kink, see [ll| . We list these coefficients up to the twentieth order in Table 
IV: 
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2.52587 x 10~ 4 
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2.66667 
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2.97161 x 10~ 5 
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1.06667 
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3.12801 x 10" 6 
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0.304762 
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2.97906 x 10~ 7 
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0.0677249 
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2.59049 x 10~ 8 
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0.0123136 
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2.072239 x 10" 9 
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1.8944 x 10" 3 
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1.5351 x 10~ 10 
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CnCPQ/a 2 "- 1 


15 


1.05869 x 10~ n 


16 


6.83027 x 10~ 13 


17 


4.13956 x 10~ 14 


18 


2.36546 x 10" 15 


19 


1.27863 x 10~ 16 


20 


6.55706 x 10~ 18 



Table IV: Seeley Coefficients 

Write now the spectral zeta functions in the form: 
Ck (s) = Ck 0; b) + B Ko (s; b) 
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(k(s) = ( K (s;b)+B K ( S ;b) 
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The incomplete Euler Gamma functions 7[z,a] are meromorphic functions of z whereas Bx (s;b) and are 
entire functions of s. The splitting point of the Mellin transform is usually taken at b = 1. We leave b as a free 
parameter for reasons to be explained later. 

Neglecting the entire parts, the zero-point energy renormalization 
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gets rid of the Cq(K) term. The contribution of c\(K) 
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is exactly canceled by the mass renormalization counter-terms: 
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We must now subtract the contribution of the zero mode: 
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Finally, the high-temperature one-loop correction to the K kink energy is: 
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In practice, truncation of the series is also necessary 
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Using formula ([21)) to calculate the one- loop kink mass shift, we admit an error of: 
Ag-A.E(6,iVo) = --3- [ y £c n (K)(T[n-l,b]+ 1 J + £ ^)r(fi-l)(H 

\n=2 ^ ' „ = at 0+ i 



r 2(n-l) 



(21) 



We offer a Figure where formula (|2~lj) has been applied for iV = 20 and several values of a. The very good precision 
of the asymptotic formula was achieved by adapting the parameter b to the value of a. For instance, we have taken 
b = 1000 for cr = 0.1, b = 100 for a = 0.3, b = 50 for a = 0.5, b = 20 for a = 0.7, & = 10 for a = 0.9, and b = 10 
for a = 1. Physically, this means that the lighter the particle mass (a 2 ) is, the longer the integration interval in the 
Mellin transform must be taken to minimize the error produced by the neglected entire parts. In practice, we have 
chosen b in each case at the frontier near the point 0a € (0, oo], where the asymptotic formula of the A'-heat trace 
departs from its exact value. 
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Figure 2: Points obtained using the asymptotic expansions for several values of a plotted on the exact curve giving the one-loop 
mass shift as a function of a. 



VII. CONCLUSIONS AND FURTHER COMMENTS 



In sum, we may draw the following conclusions: 

1. We have obtained the one- loop mass shift to the classical mass of the stable topological kink that exists in a 
massive anisotropic non-linear § 2 -sigma model. 

2. In the isotropic case, a = 1, our result agrees with the answer provided by other authors: the one- loop correction 
is twice (in modulus) the correction for the sine-Gordon kink, see Q and [l5| . 

3. Our procedure is based on the heat kernel/zeta function regularization method. The result is identical to the 
answer achieved by means of the Cahill-Comtet-Glauber formula. 

This is a remarkable fact: the CCH formula takes into account only the bound state eigenvalues and the 
thresholds to the two branches of the continuous spectrum of the Schrodinger operators that govern the field 
small fluctuations. It is essentially finite. Our computation involves infinite rcnormalizations. The criterion 
chosen to set finite renormalizations -no modification of the particle masses at the one-loop level, equivalent to 
the vanishing tadpole criterion in linear sigma models- does exactly the job. 

4. We have also derived a high-temperature approximated formula for the mass shift, relying on the heat kernel 
asymptotic expansion. We stress that we have improved a former weakness of our method. The approximation 
to the exact result was poor for light masses -non-dimensional mass < 1- in the model studied in [9(. We have 
achieved a very good approximation in this paper even for light particles by enlarging the integration interval 
of the Mellin transform and considering an optimum number of Scclcy coefficients. We believe that this is a 
general procedure, working also in models where the exact generalized zeta function is not available. 

As a final comment, we look forward to addressing the quantization procedure for: (a) Multi-solitons and breather 
modes of this model, (b) Stable topological kinks that may arise in other massive non-linear sigma models with 
different potentials, e.g., quartic, and/or different target manifolds, e.g., § 3 . 
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APPENDIX I: Kink Casimir energy and mass renormalization near the pole 



The Casimir kink energy is, see Section §. V.C: 
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but s = — A is a pole of AE c (s). To find the residue, we expand this function in the neighborhood of the pole by 
using the following results 
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where 2-F\ ' ' ' [|, 0, | , — jttt] is the derivative of the Gauss hypergeometric function with respect to the second 
argument and we made use of the fact that 2F\[\ 1 0, |, — rz— j-] = 1. 

The physical limit s = — i is also a pole of AE MR (s), see Section §. V.D: 
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(22) 



APPENDIX II: The heat kernel expansion 



Consider the Kg- and K-he&t kernels 

d 



K Ko {x,y;Q) = S{x-y) 

K K (x,y;0) = 5(x - y) 
which provide an alternative way of writing the Kg- and if-heat traces: 
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(23) 



Tr^e"^ = lim / dx K Ko (x, x; f3) 

l— too I I 



Tr L 2e- 0K = lim / dx K K (x,x; /3) 



Note that the form of the i\-hcat equation (|23|) . (j^ + K — U(x)j Kk{x,v; (3) = 0, suggests a solution based on 

the iio-heat kernel: KK{x,y; (3 — Ck (x, y; /3)Kk (x, y; (3). The density CK(x,y; (3) satisfies the infinite temperature 
condition C/f(a;,?/;0) = Inxn and the transfer equation: 



d x — y d d 2 \ 
d(3 (3 dx dx 2 J 



C K (x,y;P) = U(x)C K (x,y;(3) 



(24) 



Next we seek a power series solution, C^-(a;,y; (3) — X)^°=o c n(x,y)f3 n , of (|24|) . which becomes tantamount to the 
recurrence relations: 



0c n d 2 c n —\ 
nc n (x,y) + (x - y) — (x,y) = a " 2 (x,y) +U(x)c n - 1 (x,y) 



Ox Ox 2 v ~' ay ' -v-v-»-v~>», • ( 25 ) 

In fact, only the densities at coincident points x ~ y on the line are needed. We introduce the notation ^C n (x) = 

limz^y d Q°h y) to write the recurrence relations for the Seeley densities (and their derivatives) in the abbreviated 
form: 



{k) c n (x) = -L- 

n + k 



The (Seeley) coefficients c n (K) are the integrals over the infinite line of the densities c n (x,x), i.e., c n (K) 
J^dxcnix^). 
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